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ABSTRACT 

The MACHO project is a search for dark matter in the form of massive compact 
halo objects (MACHOs). The project has photometrically monitored tens of millions of 
stars in the Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and Galac- 
tic bulge in search of rare gravitational microlensing events caused by these otherwise 
invisible objects. In 5.7 years of observations toward the LMC some 13-17 microlensing 
events have been observed by the MACHO survey, allowing powerful statements to be 
made about the nature of the dark population in the halo of our Galaxy. A critical 
component of these statements is an accurate determination of the survey's detection 
efficiency The detection efficiency is a complicated function of temporal sampling, stel- 
lar crowding (the luminosity function), image quality, photometry, time-series analysis, 
and criteria used to select the microlensing candidates. Such a complex interdependence 
is most naturally solved using a Monte Carlo approach. Here we describe the details 
of the Monte Carlo used to calculate the efficiency presented in the MACHO 5. 7- year 
LMC results. A similar calculation was performed for MACHO's 1-year and 2-year re- 
sults. Here we correct several shortcomings of these past determinations, including (1) 
adding fainter source stars (2.5 magnitudes below our faintest detected "stars"), (2) an 
up-to-date luminosity function for the LMC, (3) better sampling of real images in both 
stellar density and observing conditions, (4) an improved scheme for adding artificial 
microlensing onto a random sample of real lightcurves, and many other improvements. 
The Monte Carlo technique presented here realistically simulates the negative effects of 
severe crowding (blending) that is a fact of microlensing surveys. 

Subject headings: dark matter — gravitational lensing — methods: data analysis - 
Stars: low-mass, brown dwarfs 
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1. INTRODUCTION 

Gravitational microlensing has become an impor- 
tant tool for the discovery and characterization of 



ing groups ( 


Afonso et al. 1999, 


Alcock et al. 1997^, 


Ansari et al. 1997 


, Crotts et al. 199E, [Jdalski et 


al. 1997 


) as there are potential lines of sight out of 



the Galactic halo (LMC, SMC, Galactic bulge, M31). 
Surveys toward the LMC have observed an excess in 
the number of microlensing events over what is ex- 
pected from known populations of stars. We have 



recently reported 13-17 microlensing events (Alcock 



et al. 1999e; hereafter A00) in 5.7 years of obser- 
vations toward the LMC and recently EROSII has 
reported two new events that they interpret as limit- 



ing the amount of halo dark matter (Lasserre et al 



2000), but are consistent with the results of A00. In 
addition, a number of candidates have been observed 



towards the SMC (Alcock et al. 1997b, 


Alcock et al. 


1999a, 


Albrow et al. 1998b, 


Palanquc-Dclabrouillc et 


al. 1998|, 


Afonso etal. 1999^ 


Afonso et al. 1998). One 



natural explanation for this excess rate is a popula- 
tion of massive compact halo objects of mass ~ 0.6M Q 
that contribute ~ 20% to the mass of our Galaxy's 
halo. However, previously unknown (or underesti- 
mated) populations of stellar lenses (e. g., in an LMC 
halo) are a distinct possibility. In order to make quan- 
titative statements about such a population an accu- 
rate determination of the survey's microlensing detec- 
tion efficiency is required. Here we describe the MA- 
CHO project's pipeline for calculating its detection 
efficiency More details on the experiment, microlens- 
ing terminology, analysis, and interpretation can be 
found in the companion paper A00, and detailed re- 
views of microlensing in general are given by Roulet 
& Mollerach 1997 and Paczyhski 1996. 

The detection probability for individual events de- 
pends on many factors, e.g., the 3 event parameters 
A max , t, t max (maximum magnification, Einstein- 
diameter crossing time and time of peak magnifica- 
tion, respectively), and the unlensed stellar magni- 
tude, as well as the observing strategy and weather 
conditions. Such a complicated dependence is most 
naturally solved using a Monte Carlo technique. We 
may simplify the dependence by averaging over the 
known distributions in A maxi t max , the stellar mag- 
nitudes, and the knagm l±me=sanipiing and ateathcr 



corJdilians, to derive efficiency a function nr| ly 
of event timescale, £(t). 



Given an efficiency £ (t) we may compute the ob- 
served optical depth as, 



7T 1 
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E 



u 



£{U) ' 



(1) 



where E is the survey's exposure in object-years, U is 
the duration of event i and £(t) is the efficiency for de- 
tecting microlensing with duration t (see also Alcock 
et al. 1996 and 1997a; hereafter A96 and A97, respec- 
tively). The optical depth is a function of the amount 
and distribution of mass along the line of sight and 
is independent of the lens masses or velocities. How- 
ever, we may extract more information if we assume a 
model of the Galactic halo. With the use of a model, 
a likelihood analysis may be performed to determine 
the most likely lens mass m and mass fraction / of the 
halo model. A model yields the distribution of event 
durations dT/di ( Alcock et al. 1995| ), which, when 



combined with the efficiency and integrated over all 
possible durations, predicts the number of events one 
expects to observe from such a Galactic halo com- 
posed entirely of MACHOs. 



N, 



E 



s(t)^di. 

dt 



(2) 



An identical method may be outlined for stellar pop- 
ulations as in A00. A likelihood estimator may be 
constructed from the observed events and the model- 
predicted number of events to determine simultane- 
ously the most likely lens mass m and mass fraction / 
of the halo. Note the vital role the efficiency plays in 
connecting the world of models to that of the survey. 

The efficiency is a strong function of the tempo- 
ral sampling of the survey. Events with very short 
durations (or very long durations) are unlikely to be 
detected as they easily fall 'in-between' observations 
(or for long durations, extend through the data win- 
dow) . A straightforward way of simulating this "sam- 
pling efficiency" while retaining realistic behavior of 
the data is to use a random sample of real lightcurves 
from the survey in which to inject artificial microlens- 
ing. By randomly generating a number of events and 
running the same time-series analysis and selection 
criteria used in the real analysis, this Monte Carlo 
can be used to determine the survey's detection effi- 



1997) 



ciency (Alcock et al. 1996, Alcock et al. 1997a, Pratt 



How does one add microlensing onto a lightcurve? 
A simple method is outlined in § 3.2 for the case where 



3 



one assumes that each lightcurve represents a single 
resolved star. In this case all of the flux represented 
by the lightcurve is magnified by the microlensing 
model A(t) (Equation ^ below). The rarity of mi- 
crolensing, however, demands crowded stellar fields 
in order to detect even a handful of events, and this 
in turn results in the blending^] of stellar light. The 
issue now becomes what fraction of the flux in the 
lightcurve (or object) is lensed? Our solution was 
first described briefly in A96 (see Pratt 1997 for a 
more thorough description) and is outlined briefly as 
follows. The method involves the use of a large set 
of artificial star tests on a sample of MACHO sur- 
vey images that are seeded with a color-magnitude 
diagram of the LMC. Each artificial star was added 
to an entire series of observing conditions and over 
15 peak magnifications. The MACHO survey's ded- 
icated photometry code SoDoPhot was then run on 
the resulting images to determine how each star be- 
haved over a range of seeing, sky, and magnification. 
The resulting photometry and photometric flags were 
stored in binary files, which we refer to as Photomet- 
ric Response DataBases (PRDBs). In this "photo- 
metric efficiency" technique each artificial star (hence 
forth referred to as a Photometric Response Func- 
tion or PRF) is used as a model for adding artificial 
microlensing onto lightcurves. In most cases only a 
fraction of the observed flux in the lightcurve is actu- 
ally lensed, and the PRF supplies us this fraction as 
a function of seeing and sky. In this way both blend- 
ing and systematic photometry effects are mimicked 
in the resulting artificial data. 

Blending introduces several problems that can af- 
fect microlensing surveys in serious ways and has been 
extensively investigated (Alcock et al. 1996, Wozniak 
19E0, [Pratt 1997| , |Han 1997j , |Han 1998| ). In addition 
to diluting the true peak magnification A max , blend- 
ing also biases the measured durations i to shorter 
values, since blended events spend less time above a 
given threshold. This biasing of t is particularly im- 
portant as the optical depth estimate is proportional 
to the average duration of the events. Another ef- 
fect of blending (which has received less attention, 
but is equally important) is that the survey's expo- 
sure in star-years is typically quite a bit larger than is 
estimated by naively counting photometered objects. 



1 Blending is when two or more star's point-spread functions 
overlap to such an extent that the photometry code can no 
longer identify individual stars, but rather recovers one photo- 
metric 'blob' or object. 



One can view this (for a fixed exposure in object- 
years) as a net increase in the efficiency which in part 
balances out a decrease in efficiency due to blending. 
All of these effects must be taken into account if one 
desires an accurate detection efficiency known at least 
to the level of, if not better than, the shot noise in- 
herent in the low number of events seen toward the 
Magellanic Clouds. 

We make several corrections and improvements to 
the Monte Carlo presented in A96 and A97. Most no- 
tably these are: (1) we add source stars to V — 24.5, 
more than 2.5 magnitudes below our faintest detected 
stars (and 2.5 magnitudes fainter than was used in 
A96 and A97), (2) we use an up-to-date luminosity 
function of the LMC, constructed with our ground- 
based photometry for bright stars (V < 19) and HST 
photometry for dim stars (V > 19), (3) we com- 
pute luminosity function normalizations separately 
for each MACHO field, (4) 10 fields of size 5' x 5' with 
widely differing stellar density were used to simulate 
photometry of artificial stars over an average of 69 
different observing conditions (only 2 fields with ~ 20 
observing conditions were used in A96 and A97), (5) 
an improved scheme for adding artificial microlensing 
onto real lightcurves, and (6) we describe a robust 
method of statistically correcting for the t bias that 
can be used to estimate the optical depth. 

In § U we briefly review the MACHO telescope, 
photometry system and time-series analysis in order 
to introduce some MACHO specific vocabulary used 
throughout the rest of the paper. § ||| outlines the 
sampling efficiency approach as a primer to the pho- 
tometric efficiency discussion in § ||. The results of 
the photometric efficiency analysis is reserved for § 0, 
and we conclude in § |^. 

2. THE MACHO TELESCOPE, PHOTOM- 
ETRY SYSTEM, AND TIME-SERIES ANAL- 
YSIS 

The MACHO project has full-time use of the 1.27m 
telescope at Mount Stromlo Observatory, Australia, 
for a period of 8 years starting July 1992. The tele- 
scope was re-commissioned especially for this project, 
and a computer-controlled pointing and drive sys- 
tem was installed. A system of corrective optics has 
been installed near the prime focus, giving a focal 
reduction to //3.9 with a 1 degree diameter field of 
view. A dichroic beam splitter and filters provide si- 
multaneous images in two passbands, a "red" band 
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(~ 590 - 780 nm) and a "blue" band (~ 450 - 590 
nm). Two large CCD cameras are employed at the 
two foci; each contains a 2 x 2 mosaic of 2048 x 2048 
pixel Loral CCD imagers. The pixel size is 15/jm, 
which corresponds to 0."63 on the sky, giving a sky 
coverage of 0.7 x 0.7 square degrees. Details of the 
camera system are given by Stubbs et al. (1993) and 
Marshall et al. (1994) while details of the telescope 
can be found in Hart et al. (1996). 

The survey's photometry code SoDoPhot is a PSF 
fitting algorithm based on Z?oP/jo£0 (3chcchtcr, Ma- 



teo | fc Saha 1993 ). It runs in two different modes. 
The first mode, a template generation mode, is de- 
signed to run on a pair of red and blue images or 
"chunks." (Each CCD is broken up into 16 roughly 
5' x 5' somewhat overlapping regions called "chunks" 
for the purpose of photometric reductions.) The red 
chunk is reduced first in a manner similar to a stan- 
dard DoPhot reduction. Next the blue chunk is re- 
duced using the results from the red chunk to warm- 
start the reduction, which improves the star matching 
between colors. This produces a set of pre-templates 
for this chunk pair. Once a set of pre-templates for 
an entire CCD (16 chunks) has been generated the 
final templates are constructed by enlarging the pre- 
templates and including stars that lie as much as ~ 90 
pixels outside the nominal chunk boundaries. This 
enlargement lessens negative effects on the routine 
photometry due to incorrect telescope pointing. The 
final set of template files contains a list of detected ob- 
jects, their positions (corrected to an airmass of one) 
and their template magnitudes. A set of template 
files need be generated only once and is created using 
images with better than average seeing and dark sky 
conditions. 

All other images are processed in "routine" mode, 
which proceeds as follows. The appropriate template 
file is used to warm-start the image by first locat- 
ing and matching 50 bright unsaturated stars. These 
stars are used to determine a seven parameter an- 
alytic fit to the PSF, a coordinate transformation, 
and a photometric zero point relative to the tem- 
plate. Then all the template stars are subtracted 
from the image using the model PSF and coordinate 
transformation; the variance estimate for each pixel 
is adjusted to allow for errors in subtraction. Next, 
photometric fitting is carried out for each star in de- 
scending order of brightness, by adding the analytic 



2 SoDoPhot = Son Of DoPhot 



model of the star back to the subtracted frame and 
fitting a two-parameter fit to the star's flux and sky 
background, with pixels weighted by inverse variance. 
The model PSF and computed position of the star are 
kept fixed. When a star is found to vary significantly 
from its template magnitude, it and its neighbors un- 
dergo a second iteration of fitting. For each star the 
estimated magnitude and error are determined, along 
with six other parameters ("flags") measuring (1) the 
object "type" (single/blended, etc.), (2) the x 2 of the 
PSF fit, (3) the "crowding," i.e., the amount of flux 
contributed from nearby stars, (4) the weighted frac- 
tions of the PSF masked due to cosmic rays, (5) the 
weighted fractions of the PSF masked due to bad pix- 
els and (6) the fit sky value. 

The set of photometric data points for each de- 
tected object are rearranged into a time series and 
combined with other relevant information (including 
seeing, average sky brightness, airmass, etc.) into 
portable binary files called sodsets. The sodsets are 
in turn passed into our time-series analysis code to 
search for variable objects and microlensing candi- 
dates (see A00 for more details of the analysis). In 
brief the analysis code calculates for each lightcurve 
a set of variability statistics, average magnitudes, er- 
ror bars, crowding, etc., and for lightcurves deemed 
interesting (level- 1 events) a five parameter fit to mi- 
crolensing is applied, where the parameters are the 
unmagnified red and blue fluxes, the peak magnifica- 
tion A maxi the time of peak magnification t maxi and 
the Einstein-diameter crossing time t. Many statistics 
describing the significance of the fit are also stored for 
the level-1 events. 

From the complete set of statistics we designed cri- 
teria that select out microlensing candidates from a 
wide background of noise induced bumps and variable 
stars. Two different sets of selection criteria (criteria 
set A and B) were used in A00, and we refer the inter- 
ested reader to this paper for full details of the two 
selection criteria sets, as we give only a qualitative 
'feel' for the criteria here. Criteria set A superficially 
resembles the criteria used in A97 and was designed 
to be somewhat restrictive in the sense that only 
good quality, high signal-to- noise (S/N) events were 
selected. The criteria are fairly tight and rely strongly 
on S/N statistics that are somewhat microlensing- 
shape specific. Only 13 candidates in the first 5.7 
years of LMC data passed criteria set A. In contrast 
criteria set B was designed to be inclusive in the sense 
that any event with a flat baseline and with one sig- 
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nificant and unique bump was included. However, 
criteria set B was somewhat vulnerable to variable 
stars that exhibit constant baselines for long periods 
of time with only one outburst, such as exhibited by 
supernovae, and possibly cataclysmic variables and 
nova (see A00 for a full discussion on how these inter- 
lopers were treated). Some 17 candidates in the first 
5.7 years of LMC data passed selection criteria set B. 

The relative looseness of criteria set B over A was 
due to a set of new statistics that allowed us to more 
accurately characterize and remove the periodic and 
quasi-peridoic variable star populations. With a ma- 
jor source of background removed the S/N level could 
be lowered. One of the main overall differences in 
the two sets of selection criteria is that criteria set A 
is less sensitive on the whole to moderately blended 
and highly blended events as compared to criteria set 
B. Also, criteria set B relies less on the microlensing- 
shape information and is less likely to be missing ex- 
otic forms of microlensing such as parallax and binary 
events. Although criteria set B has overall higher de- 
tection efficiency than criteria set A the number of de- 
tected events compensates and the results presented 
of A00 are fairly insensitive to the choice of cuts. In 
what follows we will present both sets of selection cri- 
teria and discuss their differences. 

3. THE SAMPLING EFFICIENCY 

3.1. The One Percent Database of Lightcurves 

A key element in the efficiency determination is 
the 1% database of lightcurves. This database rep- 
resents an unbiased random selection of 1% of the 
lightcurves contained in the MACHO project's LMC 
survey. It is unbiased with respect to observed mag- 
nitude, color, sampling rate, spatial distribution, and 
general data quality. The 1% database contains 
118, 645 lightcurves from the top 30 LMC fields^ and 
occupies 2.7 Gbytes of disk space. The lightcurves 
begin on day 200 (= JD - 2, 448, 623.5; 1992 July 19) 
and end on day 2267 (1998 March 17) which makes 
the data window 5.7 years long. The mean number 
of exposures per field is 719, with a range from 180 
to 1338. Approximately 13% of the lightcurves have 
information in only one passband (red or blue). Of 
these ~ 60% have only red lightcurves and ~ 40% have 
only blue lightcurves. There is a fair amount of field 



3 See A00 for a complete list of the 30 LMC fields along with RA 
and dec centers. 



overlap among the top 30 fields with about ~ 20% of 
lightcurves duplicated across fields. 

3.2. Modifying Lightcurves: Sampling Effi- 
ciency 

Here we assume that each lightcurve in the 1% 
database is a single resolved star, and that our pho- 
tometry code is "perfect" in recovering flux. We fur- 
ther assume that the measurement errors involved are 
dominated by photon shot noise. Given a lightcurve 
we proceed as follows: (1) first a robust mean mag- 
nitude m for each bandpass is computed. (2) Then 
a set of event parameters u m i„, t max , and t is gen- 
erated where, (a) u m i n is chosen as a uniform devi- 
ate from zero to the experimental threshold ut = 
1.0 (A max = 1.34), (b) t max is chosen as a uni- 
form deviate from slightly before the beginning of 
the observations to slightly after (day 190 to 2277), 
and (c) t is chosen as a uniform base 10 logarith- 
mic deviate with a duration of 1 day up to 2000 
days. (3) For each point on the lightcurve we com- 
pute the "scatter" from the mean in flux units as 
A/ = / — /, where / is a robust mean magnitude 
expressed in flux units, i.e., m = —2.5log(f). This 
scatter is assumed to be due entirely to photomet- 
ric error rather than intrinsic variation in the source 
(indeed, only ~ 1 — 2% of all LMC lightcurves show 
signs of intrinsic variability). (4) The flux / for 
this point is then "magnified" / — > A(u(t))f + A/, 
or expressed equivalently in magnitudes, m — > m — 
2Mog(A{u(t)) - 1 + lpC^- ™-)/ 2 - 5 ). (5) The flux error 
is modified as <jf — > <jf v/ A(u(t)), which translates 
into magnitude space as a — > a/ y/A(u(t)). Here t is 
the time of the observation being modified and u(t) 
is 



*(*) = 



(2(t-*o)AT 



1/2 



with 



A(u) 



u\Ju 2 + 4 



(3) 



(4) 



This is performed for all points on the lightcurve and 
for both bandpasses. We leave the remaining photo- 
metric flags (Xpsfi crowding, etc.) unchanged. 

3.3. Calculation of the Efficiency 



Using the technique of § 3.2 to modify lightcurves 



we proceed to add an event onto each lightcurve in 
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the 1% database. One complete loop through the 1% 
database we refer to as a pass. For each "event" we 
save information such as, (1) the lightcurve id number 
(field. tile. sequence), (2) a pass number, (3) a robust 
magnitude V bj and color V — i? bj of the lightcurve, 
and (4) the input event parameters A max , t max , and t. 
This bookkeeping information will later be matched 
with the output of the time-series analysis and used 
to calculate the efficiency. 

Once we have created a complete pass and stored 
the bookkeeping information we next run the time- 
series analysis on the artificial lightcurves (§ ||). The 
resulting output statistics and fits are written into a 
file and then matched on a lightcurve by lightcurve 
basis with the bookkeeping files that contain the in- 
put statistics. Hereafter all parameters labeled with 
a prime refer to measured parameters (e.g., A' max , 
t' max , and t') and un-primed parameters refer to in- 
put or true parameter values. In addition, both cri- 
teria sets A and B described in § || are applied to the 
resulting statistics, as in A00, and the results stored 
as an integer: zero = failed the criteria set, one = 
passed the criteria set. This creates a single datacube 
file containing input statistics, output statistics, and 
'cut' integers. In order to generate enough events to 
ensure adequate statistics we ran 5 passes through 
the 1% database. Though the same lightcurves were 
used 5 times each, completely different event parame- 
ters were generated for each pass. Because we stored 
a unique pass number in the datacube file we sim- 
ply concatenated all 5 datacube passes into one large 
datacube file for ease of use. 

With this datacube it is straightforward to cal- 
culate the survey's sampling efficiency. Since the 
artificial events are added with a uniform distribu- 
tion in Umin and t max , integration over these vari- 
ables is not required. Furthermore, since our as- 



sumption in § 3.2 is that each lightcurve represents 



a single resolved star, there is also no need to inte- 
grate over luminosity (since, under this assumption, 
the 1% database fairly samples the survey's luminos- 
ity distribution). We need only bin the data in t and 
simply count the number of recovered events versus 
the number of added events. The efficiency is simply 
£(t) = N rec (t)/N add (t). 



4 We ignore a small modulation (~ 5%) in the event rate due to 
the motion of the earth around the sun. 



3.4. Results of the Sampling Efficiency 

In Figure |l| we compare the sampling efficiency de- 
rived from the 5. 7- year data set for both criteria set 
A and B with the previous two data sets, A96 and 
A97. The general behavior of the efficiency with t is 
easy to understand. The gradual fall in efficiency as 
t decreases is caused by a combination of a typical 
sampling rate of 1-2 days and large gaps of 3-50 days 
that exist due to bad weather and other telescope 
down time. The sharp drop at ~ 600 days is due to 
an explicit cut requiring i < 600 days (affecting both 
criteria set). The large difference in the efficiencies 
for the two selection criteria is due to a combination 
of (1) the different values used for the A max cut (cri- 
teria set A uses A max > 1.49 while criteria set B uses 
A m ax > 1-34), (2) criteria set B uses fewer S/N statis- 
tics, most of which are looser, (3) differences in vari- 
able star cuts, and (4) various cuts on the minimum 
quality of the lightcurve (minimum number of points 
in the red and blue passbands, maximum crowding, 
etc.). 

The most striking difference between the previous 
two data sets and the 5. 7- year set is the much higher 
efficiency at long durations. Much of this difference 
is just a reflection of the longer baseline and an in- 
crease in the allowed duration (t < 600 days). How- 
ever, a significant amount of the difference lies in a 
quirk of the A97 data set. In the A97 data six of 
the densest fields^ had their lightcurves cut in half, 
roughly, due to an early generation of templates used 
to reduce the photometry for these fields. These early 
templates made use of a different star naming conven- 
tion than is currently used. Roughly halfway through 
the A97 data set the current generation of templates 
was implemented on these six fields. Since there was 
no straightforward way to cross reference the ID's of 
the stars in each different convention (indeed, the map 
was not even one-to-one), and redoing the photome- 
try was too prohibitive at the time, it was decided to 
consider as separate lightcurves the stars in these six 
fields before the upgrade and after the upgrade and 
to analysis them as such. This was handled in a self- 
consistent manner and more details can be found in 
A97. The primary effect of this was a lowered effi- 
ciency for long duration events, due primarily to (1) 
an effective cut on event peak t rnax where the fields 
were split and (2) the required minimum of 40 base- 
line points outside ± 2 1 of the peak of the event that 



5 These fields are 1, 7, 9, 77, 78, and 79 
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was used in the A97 selection criteria. The problem 
did not exist in the A96 data, thus the rather close 
behavior between A96 and A97 for long durations, 
even though the later had twice the coverage. The 
A00 data has redone all its photometry in these six 
fields with the current generation of templates and 
does not suffer from this problem. 

It is also worth noting that in the t range 1-40 days 
the A00 sampling efficiency (for both criteria) is sys- 
tematically smaller than the A96 and A97 results. In 
part this is due to differences in the selection crite- 
ria, but is also a product of including an additional 8 
less densely sampled fields into the year A00 analysis. 
Though this tends to reduce our efficiency it also in- 
creases our exposure. Also note that criteria set B is 
always above criteria set A except for durations less 
than 10 days. The slightly lower efficiency of criteria 
set B in this range is due to the larger number of re- 
quired points in the peak N pea k > 10 (criteria set A 
only required 7 or more 2a points; see A00). 

In addition to £(t) it is also interesting to investi- 
gate the efficiency as a function of other parameters, 
such as magnitude, impact parameter and stellar den- 
sity. Also of considerable interest is the t bias men- 
tioned in the introduction. However, we hold off dis- 
cussing these until after we have introduced the pho- 
tometric efficiency as this scheme more realistically 
models the systematics of blending and the survey as 
a whole. 

4. THE PHOTOMETRIC EFFICIENCY 

In the introduction some of the potential problems 
associated with the effects of blending were discussed. 
Blending causes systematic underestimation of both 
i moj and t, but also increases the number of stars to 
which the survey is sensitive. Clearly none of these ef- 
fects are taken into account by the sampling efficiency 
described in § To correct for this deficiency, and 
to be able to generate realistic microlensing events, 
we make use of a large database of artificial star pho- 
tometry as first described in A96 and A97. We have 
developed techniques that allow us to use this arti- 
ficial star photometry to inject synthetic microlens- 
ing onto 1% database lightcurves, replacing the tech- 
niques described in § 3.2. The new techniques allow us 
to empirically account for the blending problems, as 
well as other systematic photometry effects unique to 
SoDoPhot such as various correlations with seeing, 
underestimation of error bars and other systematic 



changes in the photometry flags. 

Our method involves a large set of artificial star 
tests on a sample of the MACHO survey's images 
that are seeded with a color-magnitude diagram of 
the LMC. A large set of real images that fairly sam- 
ple the survey's distributions in observing conditions 
(seeing and sky) were used to create the database of 
input photometry versus output photometry for the 
sample of artificial stars and organized into a PRDB. 
The PRDB consists of - 60, 000 individual PRFs that 
represent how a star behaves as part of a blended ob- 
ject. Each star was added to images over a range 
of observing conditions and over 15 peak magnifica- 
tions. SoDoPhot was run on the resulting images to 
determine how each star behaved over the range of 
seeing, sky, and magnification. The resulting photom- 
etry and photometric flags were stored in the PRDBs. 
Each PRF represents a rule for how a lightcurve would 
respond to the addition of flux over a wide variety of 
observing conditions. This can be represented as, 



A ICC = A lec (A, V, R, Vobj , i?obj , seeing, sky), (5) 

where A rcc is the recovered magnification, given an 
input magnification A, stellar magnitudes V and R, 
object magnitudes V bj and i? t>j, and characterized 
by an observing condition (seeing and sky). 

To implement the photometric efficiency properly 
it is important to know the underlying luminosity 
function (LF) to some limiting magnitude (in our 
case about V ~ 2 4). We discuss our determination 
of the LF in 



4.1. Next, in S 4.2 



we describe how the 
PRDBs are generated, including the images, point- 
spread functions, coo rdinate and photometric trans- 
formations. In § 4.3 we show some example PRFs 
and discuss general behaviors exhibited by the PRFs. 
Discussion of the photometric efficiency (hereafter re- 
ferred to as just the efficiency) results is left for § |[ 

4.1. The LMC Luminosity Function and Color- 
Magnitude Diagram 

It is important to know the LF of stars in the LMC 
(and to a lesser extent the color-magnitude diagram) 
in order to accurately estimate the efficiency. This is 
because the LF, along with stellar density, dictates 
how much stellar PSF overlap on the sky there is as 
a function of magnitude and thus how much blend- 
ing one expects. Also, since the survey's exposure 
is greater than the number of photometered objects 
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suggests, and this can increase the sensitivity of the 
survey to microlensing, we need a fairly accurate es- 
timate of the number of stars in our fields. 

4-1.1. The Shape of the LMC's Luminosity Function 

The MACHO fields typically go incomplete at mag- 
nitudes greater than V ~ 20 but can be complete to 

V ~ 21. Because we desire to know the LF to at least 

V ~ 24 we must resort to Hubble Space Telescope 
(HST) Wide Field Planetary Camera 2 (WFPC2) im- 
ag q data. The MACHO project haa obtained HST 



W FPC2 data with the F555W and F675W filtoro for 



three fields in the LMC bar. These are located in 
the MACHO fields 2, 11, and 13. For each field, we 
have obtained "shorts" (3-4 x 30 sec exposures) and 
"longs" (3-4 x 400-500 sec exposures), in both fil- 
ters. More details of the HST data reduction and 
analysis are contained in Alcock et al. (1999b). In 
addition we have obtained the reduced HST WFPC2 
data from 6 LMC bar fields obtained and reduced by 
Olsen 1999 (kindly provided via private communica- 
tion). These 6 fields are in the F555W and F814W 
filters and have similar "shorts" and "longs" giving 
completeness limits similar to our HST data. Olsen's 
PC was positioned on old LMC globular clusters and 
so we discarded all PC data and, in addition, must 
be careful of a small contamination from cluster stars 
that extend into the edges of the WF CCDs. Details 
of the Olsen HST reductions can be found in Olsen et 
al. 1998 and Olsen 1999. 

In Figure | is plotted the V band (F555W filter) 
LF for the 3 MACHO HST fields. Each field is nor- 
malized to have the same number of stars in the range 
17 < V < 20Q, with the fields of higher stellar density, 
field 2 and 11, normalized to field 13 (lowest density 
and lowest S/N). The shape of the 3 LFs appear con- 
sistent with one another from V ~ 17, where shot 
noise dominates, to V ~ 22.5, where differences in 
stellar density between the 3 fields cause differences 
in completeness. Figure || is similar to Figure || but 
shows Olsen's 6 fields normalized in a similar fashion 
to the lowest S/N field. Again the shapes of the 6 LFs 
appear consistent with one another between the shot 
noise on the bright end and the differing completeness 
on the dim end. 

We can quantify the above 'by eye' assessment by 
computing a \ 2 between fields along with the as- 



sociated probability of obtaining a value of x 2 /dof 
worse than the measured value. We find no signifi- 
cant difference between the 3 MACHO LFs with rel- 
ative measured x 2 /dof of 1.065, 1.253 and 1.242 for 
the combinations 2-11, 2-13, and 11-13. These cor- 
respond to probabilities of p = 0.35, 0.11, and 0.12, 
respectively, that the value x 2 /dof could be worse 
than the measured value. Of the 6 Olsen HST fields 
one stands out as anomalous in its % 2 /dof value. The 
LF of NGC 1916 is the most discrepant but appears 
to be so due to heavy differential reddening (|Qlscn ei 



1998j |Qlsen 1999| ) with relative x /dof ranging 



3 This region corresponds to where the HST fields are greater 
than 95% complete. 



from a low of 1.589 to a high of 3.777, correspond- 
ing to a range of probabilities p — 5.3 x 10~ 3 to 
1.1 x 10~ 17 . We thus discard this field from further 
analysis. The remaining 5 Olsen LFs are consistent 
with one another, with relative x 2 /dof ranging from 
a low of 1.118 to a high of 1.891, corresponding to 
a range of probabilities p — 0.26 to 1.6 x 10~ 4 . The 
slightly worse probabilities seen in the Olsen fields are 
likely due to contamination from the globular clusters 
spilling over into the WF CCDs. 

In Figure | we plot MACHO 's combined LF (3 
fields) and Olsen's combined LF (5 fields), normal- 
ized to the MACHO combined LF in the range 17 < 
V < 20. The two LFs are consistent with one an- 
other in the range 17 < V < 22 (x 2 /dof = 1.215, 
p = 0.17). Also plotted is MACHO's ground-based 
LF for field 13 (~ 350,000 objects) which is one of 
our most complete, lowest reddened fields, with ex- 
cellent template seeing and sky (we have allowed for 
a small offset of ~ 0.1 in V, see below). The shape 
of the ground-based LF in the range 17 < < 21 is 
good in comparison with the MACHO combined HST 
LF (x 2 /dof = 1.295, p = 0.10). Given the apparent 
lack of any significant difference in the shapes of the 
LFs for these 9 LMC bar fields we have chosen to 
combine 8 of them (excluding NGC 1916 because of 
high reddening) to form a LF with a fairly high S/N 
on the dim end (to V ~ 22.5) and to splice this to- 
gether with the bright end (V < 20) of the MACHO 
ground-based LF. In this way we create a high S/N 
LF that is complete in the range 15 < V < 22.5. 

There is still a question of the completeness of the 
HST LFs for magnitudes greater than V ~22.5. We 
have performed artificial stars tests on our 3 MA- 
CHO HST fields and used the resulting complete- 
ness curves to correct the MACHO LF to V ~ 24.5. 
The completeness tests probably represent a slight 
underestimate of the true LF. Because of the poten- 
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tial uncertainty in the shape of the LF for dim stars 
we have opted to create two different LFs in which 
to test for systematics. These are shown in panel 
(a) of Figure |ll| The first LF (hereafter LF X ) has 
been extended from V ~ 22 to greater magnitudes 



using a power law with a slope of 0.415 ( Alcock et al 



19!! 9b) derived from a linear regression fit between 
19.5 < V < 22. This LF probably represents an over- 
estimate of the number of faint stars. The second LF 
(hereafter LF<z) uses the completeness corrected MA- 
CHO LF as an estimate of the shape for dim stars. 
The shape of the true LF is likely to lie somewhere 
in between these two, but is probably closer to the 
completeness corrected LFi- We choose to use both 
LFs in order to estimate any systematic error induced 
by such an uncertainty. 

4.1.2. Normalization to MACHO Fields 

With the shape of the LF determined^] it only re- 
mains to find a normalization, that is, a number relat- 
ing these two LFs to the total number of stars. This 
is most convenient to do on a field by field basis as the 
stellar density changes quite rapidly across our fields 
(and also inside our fields, but we average over this). 
A normalization is calculated for each field using the 
ground-based LF for that field and the 'universal' LFs 
shown in panel (a) of Figure 0. This is illustrated in 
Figure || where four MACHO fields of widely differing 
stellar density are shown. 

First, we allow a slight offset in V magnitude that 
is field dependent and typically varies from -0.14 mag 
to 0.40 mag (Table Q). This slight offset is due to a 
combination of three effects, (1) extinction (patchy, 
even on the scale of our fields), (2) the tilt of the 
LMC's disk (maximum effect ±0.07 in our 30 fields), 
and (3) fields with poorly calibrated photometry. We 
note that the offsets are not correlated with seeing, 
sky, airmass or stellar density of the template images, 
but are strongly correlated with the color of the LMC 
clump, indicating that differences in extinction are 



7 We warn the reader that the LF in the bar appears to be 'uni- 
versal', and we use this simplifying assumption for the present 
results. The HST fields used here only represent the LF in the 
LMC's bar, while the LMC's disk LF has been shown to be 



somewhat different in that 



it. larks a 



intermediate age popu- 
lation (Geha et al. 1998, Olsen 1999). Our outer most fields 
53, 55 and 57 are farther out in the disk and could thus be in 
error, though we note that NGC 1754 is also in the disk and, 
although it shows a different star formation history than the 
bar, appare ntly its LF is still fairly similar to the bar fields 
( JOlscn 1999 ). 



the dominate cause. The offsets are derived by requir- 
ing the peaks of the clump and the tips of the giant 
branch to line up in V with the 'universal' LFs. The 
effect of extinction is to lower the number of source 
stars to which the survey is sensitive. Only 22 LMC 



fields are well-calibrated (Alcock et al. 19991) and 
these typically have the smallest offsets, while 8 LMC 
fields have only approximate calibrations and typi- 
cally have large positive offsets. Once the offsets are 
applied we next require the number stars in the 'uni- 
versal' LFs in the range 17.5 < V < 18.50 to match 
the number of objects in the ground-based LF in the 
same range. The two shifts typically align the 'uni- 
versal' LF to the observed LF quite well, both inside 
and outside the calibration bin, as seen in Figure ^[ 

Figure || displays the normalizations for each of the 
30 MACHO fields derived in the manner described 
above. Our choice of units for the normalization is the 
ratio of the number of stars to the number of objects 
brighter than magnitude V. In this case V = 22 is 
convenient as there are few objects dimmer than this 
limit. The SfO(V < 22) ratios are plotted versus the 
average density (in objects/ □') for each field as solid 
circles and labeled with the field number. The nor- 
malizations are also supplied in Table [|. Note the gen- 
eral trend of increasing S/0(V < 22) with increasing 
object density, which is simply a reflection of the fact 
that our fields become less complete in more crowded 
fields. The substantial amount of scatter seen at any 
given density is a reflection of other factors affecting 
completeness. For example, SfO(V < 22) is strongly 
correlated with template sky for any given density, 
in the sense that fields with high sky tend to be 
less complete (have higher S/0(V < 22)). Though 
S/0(V < 22) appears not to be significantly corre- 
lated with template seeing this is not surprising as 
the range in seeing for our templates is not nearly as 
large as the range in sky. 

Using our 3 HST fields we may make a direct 
determination of S/0(V < 22) by simply counting 
the MACHO objects that lie inside the WF CCDs. 
These are shown in Figure ^ as solid triangles. We 
have connected points from the same fields with a 
solid line. Note that the density shown for the field 
SfO(V < 22) normalizations is an average over the 
whole field, while for the HST S/O {V < 22) normal- 
izations it is only over the WF CCDs. In all three 



'This region corresponds to where the MACHO ground-based 
fields are greater than 95% complete. 
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cases the HST frames are in higher than average den- 
sity regions of the MACHO fields, resulting in higher 
S/0(V < 22) ratios. The overall agreement between 
S/0(V < 22) as determined using the individual field 
normalizations and as determined directly with HST 
is reassuring. 

We have also created three synthetic images as an 
additional check on the normalizations and to look 
for possible biases due to blending. Briefly, the syn- 
thetic images are 5' x 5' in size and were generated 
using empirical PSFs (see § |4.2| ) derived from the tem- 
plate images. A LF similar to LF\ above was used 
to add ~ 500, 000 stars down to V — 24 uniformly 
over the image. A uniform sky with Poisson noise 
was added to match the ADU/pixel distribution and 
then SoDoPhot was run in template generation mode 
on these synthetic images. The resulting number of 
recovered objects and sky were compared with the 
real image values simulated. Several iterations over 
the number of added stars versus added uniform sky 
were required to match the observed number of ob- 
jects and sky level. The S/0(V < 22) ratio for each 
of these synthetic images is shown as an open trian- 
gle in Figure || and is connected by a solid line to 
its corresponding field average. The densities of the 
three synthetic images are in general different from 
the corresponding field averages, based on where the 
chunk lies in the field. Note the very similar behavior 
of the synthetically derived S/0(V < 22) compared 
with the field average S/0(V < 22) values and the 
HST derived S/0(V < 22). In all three cases the 
higher density had the higher S/OiV < 22) ratio. 
Also note the much larger increase in S/0(V < 22) 
for the higher densities. This suggests a saturation 
point is reached in our densest fields whereby the ad- 
dition of more stars results in a dwindling increase in 
the number of detected objects as the stars pile on 
top of one another in the image. We see a similar 
behavior in the synthetic images where to recover a 
given number of objects one must add a proportion- 
ally larger number of stars. 

Using the synthetic images we may also estimate 
the size of a possible bias in the normalizations due 
to blending. It would be desirable for our normal- 
ization bin (17.5 < V < 18.5) to contain an equal 
number of photometered objects and real stars: i.e., 
iVobject 8 (17.5 < V < 18.5)/A stars (17.5 < V < 18.5) = 
1.0. Our synthetic images have values of 1.04±0.08, 
1.04±0.05, and 1.18±0.04, where the errors are la Pois- 
son errors. They all appear to have slightly more pho- 



tometered objects than real stars in the normalization 
bin, though two are entirely consistent with 1.0. Only 
the highest density image has a value that is signifi- 
cantly larger than 1.0 (but by 4.5<r). This has the po- 
tential of biasing our derived values of SfO(V < 22) 
to larger values. However, of the three synthetic im- 
ages two have extreme densities (as these were chosen 
to look for this effect) and are not representative of 
the average density in our fields (see Figure [jj). Be- 
cause we expect the biases to be strongest in our most 
crowded fields (and to a lesser extent worse seeing 
template images) the overall bias is certainly much 
less than the 18% seen in the extreme artificial chunk, 
and probably even less than the straight average of 
~9%. Another estimate^ puts the possible bias at 6- 
8%. Given the fairly small size of the effect (< 10%) 
we make no correction for this bias, but include it in 
our error budget (see § 5.3). 

The errors in S/ O (Table 0) were estimated as fol- 
lows: (1) we allowed an uncertainty in the offsets 
of 0.1 mag, based on the scatter in the offsets be- 
tween fields, and (2) based on the estimates of the 
maximum blending bias discussed above, the normal- 
izations could be at most 10% too high and prob- 
ably only 5% too low. We then propagated these 
two uncertainties through the normalization proce- 
dure to produce the final errors in the table. The 
average error in S/O per field is ~ 20%. The den- 
sity (objects/ □') weighted mean number of stars to 
objects with V < 24 is S /Q{V < 24) = 17.05 ±3.0 
for LFi and S/0(V < 24) = 10.83±1.9 for LF 2 . The 
normalization of each field amounts to correcting the 
MACHO nominal exposure in object-years to an ex- 
posure in star-years. The uncertainty in the normal- 
ization is moderately large and translates directly into 
a moderately large uncertainty in the efficiency as dis- 
cussed in § H 

Jf.,1.3. LMC Color-Magnitude Diagram 

Since we simultaneously image in both a red and a 
blue passband we require a color-magnitude diagram 
(CMD) from which to draw our artificial stars. Such 
a CMD was created for precisely this purpose and 
its construction is described in detail in Alcock et al. 
(1999b). Briefly, we spliced the bright end (V < 18.7) 



9 This second estimate uses a simple physical model of blending 
that fits a model through the three values of iV bj C cts (17.5 < 
V < 18.5) /jV stars (17.5 < V < 18.5) versus 7V stors (17.5 < V < 
18. 5) (seeing) 2 and extrapolates the results to all 30 fields. A 
net bias of 6-8% is estimated. 
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of a CMD derived from — 9 million objects (Alcock ei 



al. 1999d) from our 22 calibrated ground-based fields 
onto the dim end (V > 18.7) of a CMD derived from 
our combined 3 HST fields (which is scaled by relative 
sky area to the 9M CMD). The HST CMD was cor- 
rected on the dim end (V > 22) using a completeness 
function that corrected the data to a slope of 0.415 
in log(dN/dV) versus V. A small amount of editing 
was also required to remove isolated high pixels (due 
to the scaling up of the lower S/N HST data) and 
bright foreground stars. The CMD used here need 
not be accurate (though we endeavored to make it 
so) as its only purpose is to draw artificial stars from 
a realistic distribution in color. We expect no depen- 
dence of the efficiency on color (and see none) and so 
henceforth discard any color dependency. 

4.2. The Photometric Response DataBase 

To create the photometric response database we 
selected 10 sub-images, approximately 5' x 5' in size 
("chunks"), that span the range of stellar densities 
observed in the 30 MACHO fields (~ 93 to - 317 
objects/D'; see Table [l]) and that lie at a range of 
distances from the optical center of the bar (from 10 
arcmins to 3.2 degrees). For each of these 10 chunks 
we extracted from the MACHO archive a range of ob- 
serving conditions that fairly sample the survey's dis- 
tribution in seeing and sky. Table || lists some relevant 
parameters for each chunk, including MACHO field 
and chunk id, template observation number, density 
in objects/D', number of observing conditions N conc i, 
number of stars per grid Nprf /grid (see below), and 
the total number of useful PRFs. The mean number 
of observing conditions for these 10 chunks is 69, a 
factor of 3 more than used in A96 and A97. 

In order to add artificial stars to these images we 
require knowledge of the point-spread function (PSF) 
for each image. We also need a photometric and coor- 
dinate transformation (CTR) that maps the photom- 
etry and position of the artificial stars to some refer- 
ence image, in this case the template images for each 
chunk. We generated empirical PSFs using a PERL 
script that automates the astronomical photometry 
package DaoPhot/AllStar ( [Stetson 1992] ). The script 
iterates multiple times with a series of simple ana- 
lytic PSFs searching for stars, fitting, and subtract- 
ing neighbors to candidate PSF stars on each iteration 
while using a progressively more realistic model of the 
PSF. The final PSF is a Lorentz-Gaussian analytic 
model with an empirical look-up table of corrections. 



PSFs for each of the ~ 1600 chunks are generated in 
this fashion, and are of higher quality than the simple 
modified Gaussian analytic model used by SoDoPhot. 
We visually inspected a number of PSF subtracted 
images and were quite satisfied with the subtraction. 
AllStar photometry, using these PSFs, was run on 
all ~ 1600 chunks and the results fed into a custom 
code that searches for and computes both photomet- 
ric and coordinate transformations. The code is based 
on the Groth (1986) algorithm and incorporates a 
color dependent term in both photometric and co- 
ordinate transformations. The inclusion of color in 
the transformations allows for the effects of airmass, 
differential refraction, and a known systematic pier 
side CCD effect in the MACHO data ("blue jitter"; 
see Alcock et al. (1999b)). From the residuals to the 
fits we find the coordinate transformations are good 
to — 0.18 pixels and the photometric transformations 
to -0.08 mag. 

Grids of artificial stars are added to each of the 
10 chunks over all — 69 observing conditions using 
the appropriate PSFs and CTRs for each observing 
condition. The grids contain from 64 to 156 artificial 
stars (Table |J) positioned on a pseudo-random spa- 
tial grid such that stars are never closer than — 21 
arcsecs (or approximately — 7 FWHM in median see- 
ing) and never closer than ~ 20 arcsecs to a bound- 
ary of the chunk. This is to avoid missing too much 
data because of telescope pointing errors. The mag- 
nitudes of the artificial stars, V and R, are drawn 
from the CMD of § |4.1.3| in the range 16.5 <V < 24.5 
and -0.5 < V - R < 1.5. The CMD is sampled in a 
square-root fashion to ensure uniform statistics over 
the large range in luminosity. Each star in the grid 
is added over a range of 15 peak magnifications A = 
0.8, 1.0, 1.1, 1.2, 1.34, 1.5, 2, 3, 4, 5, 10, 20, 30, 40, 50. In 
each image the peak magnifications for each star are 
staggered so as not to add too many highly magni- 
fied stars to a given image. To cover as much of each 
chunk's image plane as possible, 60 different grids are 
created for each chunk. This corresponds to cover- 
ing ~ 30 — 70% of the image plane (depending on the 
chunk) and represents a statistically significant num- 
ber of possible photometric conditions in which a star 
could reside. 

All stars in a grid are first added with a magni- 
fication of one to the template image of the chunk 
and SoDoPhot is run in template generation mode 
on these template images to create the template files. 
With the template files created, we next loop through 
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all observing conditions and all possible peak magni- 
fications, adding the grid to the corresponding ob- 
serving condition using the appropriate PSFs, CTRs, 
CCD gains, and with a staggered magnification for 
each star. SoDoPhot in routine mode is run on 
each image created and the results are organized and 
stored in PRDB files. We chose to store all photom- 
etry on the three spatially nearest recovered objects 
to each artificial star in order to map out where the 
added flux goes. For each grid of stars a PRDB file 
with 17.4 Mbytes of photometry is created. In to- 
tal the PRDB contains 10.4 Gbytes of photometry on 
196,740 recovered objects. 

The total number of artificial stars added to all 
10 chunks over the 60 grids per chunk is 65,580. Of 
these 54,981 artificial stars contain enough reduced 
data (greater than 50% in each passbands) to be use- 
ful as PRFs. This is a factor of five more PRFs than 
was used in A96 and A97. Approximately 16% of the 
PRFs were lost due to a combination of (1) falling 
inside an obliterated region caused by a saturated 
star (~ 2%), (2) missing data entirely in either the 
red (~ 3%) or blue (~ 6%) bandpass, and (3) falling 
near an obliterated region, CCD defect, or chunk edge 
such that on average more than 50% of the photome- 
try was missing (~5%). Since we are only concerned 
with how SoDoPhot responds to added flux the loss of 
poorly determined PRFs is not a concern. The num- 
ber of single bandpass lightcurves, CCD defects, and 
other missing data are preserved in the 1% database 
of lightcurves. 

4.3. The Photometric Response Functions 

Initially we were concerned with how SoDoPhot 
divides up flux between two or more closely spaced 
stars, and which neighbor is most 'sensitive' to the 
added flux. For example, in A00 two lightcurves 
(events 7a and 7b) are, in fact, the same event. 
Event 7 was bright enough and in a locally crowded 
enough region that some of the flux from the pri- 
mary (7a) bled into a secondary neighbor (7b) caus- 
ing a spurious detection. Although event 7b was re- 
moved from the list of microlensing events and had 
no adverse effect on the results of A00, it underscores 
the need to investigate what effect multiple crowded 
neighbors have on our efficiency. To investigate this 
we stored the photometry for the three spatially near- 
est objects to each artificial star in the PRDBs. We 
found that in the vast majority of cases (~ 97%) the 
nearest neighbor snatched most of the added flux and 



was the most 'sensitive'. This is reassuring, both for 
simplicity's sake and because our previous work (A96 
and A97) had implicitly assumed this behavior. We 
illustrate this as follows. 

For each neighbor we define a recovered magnifica- 
tion, A rec (iJ) = f(i,j)/f, where f(i,j) is the flux of 
the neighbor for observing condition i and peak mag- 
nification j. Here / is the baseline flux for the neigh- 
bor (that is the flux averaged over all i with j = 
corresponding to A = 1.0). A scatter plot of A rec ver- 
sus A for all i shows an excellent linear relationship as 
expected (see below), with a fit slope m = dA rec /dA 
that is generally between and 1. This slope corre- 
sponds to the average blend fraction. We constructed 
a measure of sensitivity to added flux and computed 
this measure for each of the three neighbors to each 
artificial star. The measure we used was the fit slope 
m divided by a mean relative error bar a for the base- 
line, i. e., measure = m/er. This measure ensured that 
neighbors with large slopes and small error bars are 
counted as the most sensitive, while in the case of 
two neighbors with equal slopes the neighbor with 
the smallest error bars is counted as the most sensi- 
tive. This is likely to be the case as both selection 
criteria set A and B make use of a signal-to-noise cut 
requiring the magnification to exceed some multiple 
of the mean relative error bar (see § || and A00) . 

For each artificial star, the three nearest neigh- 
bors are ranked according to the measure, and a cut 
A m ax > 1-34 was applied to each neighbor to ensure 
that the PRF was not 'junk' (here A max is the maxi- 
mum magnification of the recovered neighbor). Only 
in ^3% of the cases was either the 2nd or 3rd nearest 
neighbor more 'sensitive' than the nearest neighbor. 
We visually inspected a large number of these cases. 
In most cases there were two neighbors about equally 
spaced from the artificial star with neighbor 1 (the 
closest) dimmer than neighbor 2. Neighbor 2 was 
more 'sensitive' for two reasons, (1) being brighter it 
had the smaller relative error bar and (2) due to a 
slight systematic bias in SoDoPhot the brighter stars, 
which are reduced first, can pirate flux from the wings 
of nearby dim stars. Since only ~ 3% of cases are in 
error due to this effect, we restrict our attention to 
only the nearest neighbor. 

A scatter plot of A rec versus A over all observ- 
ing conditions is shown for four PRFs in Figure 0. 
Note the excellent linear relationship fit by the solid 
line. The dashed line is the simple analyt ic re sponse 
function used in the sampling efficiency (§ 3.2) and is 



13 



plotted here for comparison. The panels correspond 
to (a) a unblended PRF (to = 1.0), (b) a somewhat 
blended PRF (to = 0.77), (c) a moderately blended 
PRF (to = 0.65), and (d) a heavily blended PRF 
(to = 0.35). The scatter at each A is composed of 
~ 69 observing conditions and is typically well corre- 
lated with seeing in the sense that worse seeing in- 
duces slightly larger magnifications. This is easy to 
understand, since worse seeing implies the PSFs over- 
lap more and thus more flux can be contributed to 
the nearest neighbor by a magnified star. A few cases 
of very high positive correlation and even negative 
correlation exist and correspond to cases where (1) 
the artificial star is not directly detected and lies a 
moderate distance from the nearest neighbor so that 
when it is magnified it only affects its nearest neigh- 
bor in poor seeing and (2) the artificial star lies almost 
equidistant between the two nearest neighbors and in 
good seeing contributes to the nearest neighbor but 
in bad seeing contributes more to the second nearest 
neighbor. 

To investigate how the relative error bars a (error 
bars expressed in magnitudes) behave versus A we de- 
fine an effective error 'de-magnification' as A a (i,j) = 
cr(i,j)/<j(i,j — 0), where u{i,j) is a for observing 
condition i and peak magnification j. We normalize 
the error bar magnification separately for each ob- 
serving condition i because a is in general highly cor- 
related with seeing (becoming larger in worse seeing) . 
Notice that A a (i,j) — ► 1.0 for all i when j ^ 
(that is A —>■ 1.0). Figure || is a scatter plot of 
A a versus A for the same set of PRFs as shown in 
Figure ff| Over-plotted as a solid line on each panel 
is the purely Poisson behavior used in the sampling 
efficiency (A a = 1.0/ \/[A)). The behaviors of the 
PRF's relative error bars show little resemblance to 
this purely Poisson approximation. Rather the PRFs 
behave more like A a = 1.0/A (dashed line in the fig- 
ure), which can be understood as follows. 

In our most common case the noise is dominated 
by a combination of the Poisson noise in the sky plus 
extra Poisson noise added due to neighboring stars 
which have been subtracted during the photometry 
reductions. These are both independent of magnifica- 
tion, so the error in flux should not depend on magni- 
fication. That is, the error in linear flux units should 
not change at all, while the relative errors should be- 
have as ~ 1.0/ flux or, as in our case ~ 1.0/A This is 
the pattern we see in our PRFs. However, very bright 
stars typically fall into the purely Poisson limit. Heav- 



ily blended PRFs can also approach the Poisson limit, 
and in some cases their relative error bars change very 
little as the magnified artificial star is only a small 
perturbation on the much brighter neighbor. 

The broad range in behavior of the PRFs is illus- 
trated in Figure ^, which is a scatter plot of fit slope 
Tn dA rGC J dA versus magnitude V for ~ 59, 000 
PRFs. The amount and distribution of blending is 
immediately apparent in this plot, with the artificial 
stars being strongly bifurcated towards either being 
recovered within ± 10% of their input flux or being 
blended by greater than ~ 90%. However, there are 
still a substantial number that are recovered at in- 
termediate blend fractions, and as the artificial stars 
shown here represents a square-root sampling of the 
CMD (thus under- weighting dim stars; see § 4.2 ) 



the corresponding results for a linear sampling would 
show proportionally even more blending. The left or- 
dinate is the fit slope to and the right ordinate is 
scaled to the maximum recovered magnification that 
each PRF contains (due to our maximum input mag- 
nification of 50). For example, an artificial star with 
to = 0.5 would contain a maximum recovered mag- 
nification of A rec = 25 in the PRF (with any higher 
magnificati ons n eeding to be extrapolated beyond this 
point, see § 4.4). 

Also plotted in Figure ^ is a family of five smooth 
curves that correspond to model m(V). Each curve in 
the family assumes a certain amount of blended flux 
corresponding to Vbiend — ^8, 19, 20, 21, and 22 mags 
in the figure. For example, a star blended with 20 
mags of additional flux (the center curve) would fall 
at to = 0.7 if it was V = 19, m = 0.5 if it was V = 20, 
and to = 0.25 if it was V — 21. Note how the family of 
curves brackets the scattered points fairly well. The 
solid horizontal line labeled A — 1.75 illustrates that 
stars fainter than V ~24 rarely, if ever, are magnified 
greater than this limit. As a consequence, if one de- 
sired to add stars fainter than V ~ 24, one would also 
need to add them at magnifications greater than 50. 

The PRFs also allow us to model how the photom- 
etry flags (crowding, Xpsf' n * sky etc.) are handled 
by SoDoPhot under various observing conditions and 
peak magnifications. We briefly summarize the note- 
worthy effects here. The crowding parameter, which 
is a measure of the amount of contaminating flux from 
nearby neighbors, not surprisingly, is highly corre- 
lated with seeing. As such we define a multiplicative 
'magnification' in a fashion identical to the relative er- 
ror bars discussed above. In most cases (>80%) there 
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is little or no variation of crowding with A. How- 
ever, in a few cases the crowding parameter increases 
smoothly with A. The Xpsf * s a -^ so strongly corre- 
lated with seeing, but in the sense that the fit to the 
PSF is worse in better seeing. This anti-correlation 
is mostly due to the photometered object being made 
up of multiple overlapping PSFs, which are smoothed 
out in bad seeing. Unlike the crowding parameter, the 
Xpsf 1S anti-correlated with A, in the sense that high 
magnifications result in poor fits to the PSF. The fit 
sky parameter is not significantly correlated with ei- 
ther seeing or sky, and in only a few cases (< 10%) is 
the fit sky well correlated with A, usually in the sense 
that it is higher for larger peak magnifications, but 
not always. Neither the missing pixel or cosmic ray 
flags significantly correlate with seeing, sky, or mag- 
nification. The object type is more complicated, but 
fortunately for our purposes is not important. How 
the object type changes with seeing and sky is al- 
ready properly handled in the 1% database and how 
it might change with magnification is inconsequential 
as long as it remains a valid type. 

4.4. Modifying Lightcurves: Photometric Ef- 
ficiency 

We now replace the sampling efficiency rules (§ |3.2| ) 
for modifying 1% database lightcurves with a new set 
of rules derived from the PRFs. The PRFs offer us 
an empirical set of rules that realistically incorporate 
blending and the many systematic photometry effects 
observed in SoDoPhot that were described in the last 
section. 

To add microlensing onto a lightcurve with the 
PRFs we must first match a lightcurve to a PRF. This 
is most consistently performed using only the avail- 
able measured parameters of the lightcurves, such 
as magnitude, color, average crowding, average er- 
ror, etc. We have limited ourselves to three parame- 
ters, Vobj, V — Robj, and crowding for the match. A 
match in average error was not chosen as it is already 
highly correlated with magnitude. Instead we chose 
to match in crowding for the following reasons, (1) the 
average error is also correlated with crowding (inde- 
pendent of magnitude), (2) the crowding parameter is 
a natural measure for parameterizing the image plane 
(i. e., regions of high stellar density versus low stellar 
density) as there is a strong correlation (r = 0.95) of 
the average crowding in MACHO fields with stellar 
density, and (3) of all the measured parameters, we 
believe crowding to be the most likely to be connected 



with blending. Though we see no strong correlations 
between blending and crowding, we do see a bifurca- 
tion reminiscent of Figure ^ with a slight preference 
for strongly blended events to also be highly crowded. 

To match our two large databases (1% database 
and PRDB database) as uniformly as possible we 
use the following scheme. First, the 1% database 
of lightcurves is sampled uniformly because of the 
importance of correctly weighting the temporal sam- 
pling and the unbiased nature of the database. For 
each lightcurve in the 1% database a robust mean 
magnitude V 00 j, color V — R 00 j, and crowding C bj 
are computed. The PRFs are binned in three di- 
mensions corresponding to their baseline magnitude 
Vphf, color V — Rprf, and crowding Cppp, with 
bins sizes of ~ 0.2 magnitude (Nuns = 28), ~ 0.1 
color (N bi ns = 17), and ~ 30 crowding (N bins = 
7). The bin sizes are constant over the most dense 
regions of their respective distributions, but grow 
slightly in size near the edges to accommodate for 
sparseness. The mean number of PRFs in each bin 
is ~ 17. However, due to the shape of the CMD 
~ 51% of these bins are empty (as they should be) 
and thus the mean number of PRFs in occupied bins 
is ~ 34, with the typical occupied bin containing 12- 
20 PRFs. With the large number of PRFs employed 
it is unlikely (~ 1%) that a lightcurve will encounter 
an empty bin. When this does occur we simply match 
to the nearest non-empty bin (while holding crowd- 
ing constant) which is never more than one bin away. 
The PRF bin that best matches the lightcurve is then 
selected and a PRF is randomly chosen from this bin. 
In this way we sample the PRFs as uniformly as pos- 
sible. 

Once a PRF has been matched to a lightcurve, 
event parameters A max , t max , and t are generated as 
in § 3.2, The event is now added onto the lightcurve 



using the rules from the PRF. Each observation on the 
lightcurve with A(t) > 1.005 is matched to the clos- 
est corresponding observing condition, i, in the PRF. 
In practice this amounts to minimizing the quantity 
As = (Aseeing) + (Asky) 2 to determine i, where 
Aseeing is the difference between the observation's 
seeing and the seeing in one of the ~ 69 observing 
conditions contained in the PRF, normalized to the 
maximum range of seeing. Asky is similarly defined 
and normalized. 



10 The crowding parameter ranges from to 255 and is logarith- 
mic in nature. 
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Having determined the observing conditions i, a 
two point linear interpolation is used to compute A rec 
based on the low point Ai a — A rec (i,j) and high 
point Ahi = A rec (i,j + 1) that bracket A(t) (or the 
lowest two points for an extrapolation). To preserve 
as much of the intrinsic lightcurve scatter as possible a 
technique similar to that outlined in § |3.2| is employed 
in magnifying the flux. That is each magnitude is 
modified asm -> m-2.5log(A rec (u(t)) - 1 + i()(™- m )/ 2 - 5 ) 
The relative error bars are modified as a — > A a a, 
where er is the relative error bar for this observation 
and A a is also computed via a linear interpolation 
between A i Q a nd Am from the tabulated values of 
A<r(i,j) (§ |T3|). This technique preserves the intrin- 
sic size of the lightcurve's relative error bars and re- 
alistically modifies how they respond with A over the 
various observing conditions. 

The PRFs also allow us to modify the photom- 
etry flags recorded by SoDoPhot. As discussed in 
§ O some of these show clear signs of systematic 
behavior with magnification. We chose to modify (1) 
the crowding parameter, which represents the amount 
of flux contributed from nearby stars, (2) the Xpsf? 
which tends toward worse fits in better seeing and 
high magnification, and (3) the fit sky value. We do 
not modify (4) the object type flag, (5) the weighted 
fractions of the PSF masked due to cosmic rays, and 
(6) the weighted fractions of the PSF masked due to 
bad pixels. 

The flags are modified as Flag — > Flag — dFlag, 
where Flag is the respective flag (crowding, Xpsf> or 
fit sky) and dFlag is a linearly interpolated difference 
between the PRF's flag value at A and at A = 1.0, 
computed separately for each observing condition i. 
Because the flags are logarithmic in nature the mag- 
nification translates into an additive term, much like 
in the case of magnitudes. This technique preserves 
the lightcurve's intrinsic flag values while realistically 
altering them as a function of A and observing con- 
ditions. 

A sample lightcurve modified using the four PRFs 
displayed in Figures |^ and || is illustrated in Figure [!(]. 
In each panel the input event parameters were the 
same, A max = 2.0 and i = 80.0 days. The pan- 
els correspond to blend fractions (a) m = 1.0, (b) 
to = 0.77, (c) to = 0.65, and (d) to = 0.35. Note 
the dramatic drop in magnification A max and rel- 
ative shortening of the duration i as the event be- 
comes more blended. Also note the chromatic behav- 
ior exhibited in PRF (c) due to the difference in color 



between the lensed star and the blended flux. We 
have visually inspected a large number of artificial 
lightcurves and compared the modified portions with 
real events seen towards both the LMC and Galactic 
Bulge. The comparison of magnification, error bars, 
and flags and their correlations with seeing and mag- 
nification are all quite satisfactory. 

5. RESULTS OF THE PHOTOMETRIC EF- 
FICIENCY 



A procedure similar to that outlined in § 3.3 is 



used to compute the photometric efficiency, with one 
important difference. In § 3.3 the distribution in lu- 



minosity of the stars was assumed to be the same 
as that observed by the survey (i.e., uncorrected for 
incompleteness) since all photometered objects were 
counted as resolved stars. Since we are now adding 
events onto lightcurves using an underlying luminos- 
ity distribution we need to integrate over this distribu- 
tion, as the observed distribution is clearly incorrect 
(§ 4.1). In principle, one could quantify the efficiency 



as a function of both duration t and luminosity V, 
and use this directly. However, in practice we cannot 
unambiguously determine the unlensed luminosity (at 
least without additional intensive follow-up photome- 
try for blending fits or photometry from space). Since 
the LMC's LF is well known to V ^24 we opt to in- 
tegrate out this variable. 

To generate adequate statistics we make 10 passes 
through the 1% database of lightcurves, matching 
each lightcurve to a P RF, which inturn is used to 
modify the lightcurve (§ 4.4). Bookkeeping informa- 
tion (lightcurve id, PRF id, V, V - R, V ohh V - R ohh 
A m ax, t max , t, etc.) for all lightcurves is stored and 
matched to the corresponding statistics (A' max , t' m , 
i', etc.) generated by the time-series analysis (§ |J). 
Selection criteria set A and B are applied to these 
statistics and the results stored as integers in the dat- 
acube. 

5.1. Efficiency 

The datacube is binned in two dimensions, V and 
t. Recall that events are added uniformly in t max and 
u mini and so these are averaged over. We choose 100 
bins in V of size 0.1 mag in the range 15 < V < 25, 
and 24 bins in t that are logarithmically spaced in the 
range 1.0 < i < 2000.0 days (Table |). Our artificial 
events are added uniformly in log(t), so bins loga- 
rithmically spaced insures equal numbers of events in 
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each t bin. For each bin an efficiency is computed: 
e(t,V\= N rec (i,V)/N add (t,V) where N rRC (t,V) and 
N a dd{i, V) are the number of events that pass the se- 
lection criteria set and the number of added events 
in the bin of t and V, respectively. It is worth men- 
tioning that because e(i, V) is computed separately 
in each V bin, this function is independent of the 
CMD used to seed the artificial star tests. From 10 
passes through the 1% database, ~1.2 million artifi- 
cial events are generated and of these only ~ 43, 000 
pass criteria set A and ~ 52, 000 pass criteria set B. 
The mean number of added events per bin is 500 and 
the mean number of recovered events per bin is 18 
(criteria set A). 

The function e(t, V) for criteria set B is shown in 
Figure [ll] as a contour plot in panel (b) . The contours 
correspond to efficiencies of 0.001, 0.01, 0.1, 0.2, 0.3, 
0.4, and 0.5. The gross behavior with t and V is ap- 
parent, with a broad peak in efficiency over the ranges 
30 < i < 200 days and 17 < V < 19.5. The function 
e(t, V) falls off rapidly for magnitudes smaller than 
V ~ 17 due to an explicit cut V > 17 (criteria set 
B), but the more gradual fall off for dim stars is a 
natural consequence of fainter, lower S/N events. The 
gradual drop in efficiency for short duration events is 
caused by the sampling issues discussed in § 3.4. The 



sharp cut at long durations (t < 600 days) seen in the 
sampling efficiency has here disappeared, although a 
remnant can still be seen for bright events. The cause 
of this, of course, is blending. Bright events are far 
less likely to be heavily blended, and thus are un- 
likely to be recovered with t's longer than 600 days. 
Events that are fainter are more commonly blended 
(fit t's much shorter) and can be scattered below the 
cut at 600 days. A corollary to this is that intrinsi- 
cally short events are unlikely to be detected on faint 
stars, which are typically blended, as is also seen in 
Figure [ll], panel (b). 

Given the uncertain knowledge of the source star's 
luminosity, the somewhat noisy function e(i, V), and 
the fairly well constrained shape of the LMC's LF, 
we have opted to integrate out the variable V in the 
function e(t,V). Although the shape of the LMC's 
LF is well known, the overall normalization to our 
fields, that is the numb er o f stars with V < 24 in our 
fields, is less certain (§ 4.1). This uncertainty in nor- 
malization translates into an uncertainty in the sur- 
vey's exposure in star-years. However, it is expected 
that the efficiency times the exposure will converge at 
some magnitude Kstop- Furthermore it is convenient 



to refer to our 'exposure' in units of objects-years (the 
number of lightcurves monitored times the duration of 
the survey) as this number is well known. We opt to 
move the (uncertain) normalization, that translates 
the 'exposure' in object-years into exposure in star- 
years, into the efficiency. We can understand this as 
follows. Assume a general distribution of event dura- 
tions dT/di for some Galactic model. The number of 
expected events N exp is just, 



N exp — E* 



o at 



(0) 



where the exposure here, £7*, is in star- years and the 
efficiency is calculated from e(t, V) as, 



£.(*) = 



e(t,V) ct>{V)dV. 



(7) 



Here V s t op is some cut-off magnitude where the inte- 
gration is stopped and 4>(V) is the LMC's LF normal- 
ized such that, 



<j){V)dV = 1.0. 



(8) 



The exposure, E* , in Equation g must be related to 
the observed number of objects, A Q bj, as, 



E* 



[S/0(V stop )]N obi T, 



(9) 



where T is the time-span of the survey and S/0(V s t op ) 
is a scaling factor that converts the average num- 
ber of objects seen in the survey to the actual num- 
ber of stars, down to the cut-off magnitude V^op- 
S/O(y stop ) may be estimated from, 



j-v.t= P dy 
J VBt ° P $o hi (V ohi )dV ohi 



S/0(V stop ) 



(10) 



where $ bj (K>bj) is the ground-based luminosity func- 
tion and is normalized to the total number of ob- 
jects observed by the survey: i.e., $ bj(K)bj) = 
A r obj</ , obj(V^bj)- Similarly <j>( V) is the true under- 



lying LF normalized as in § 4.1.2| . The normal- 
ization S/0(V atop ) was estimated for each field in 
§ [O] (Table [l] and Figure ||) for two cut-off magni- 
tudes V^top = 22 and V^top = 24 and for two LF 
(LFi and Li 7 ^)- We make a simple re-definition; let 
£{i) = [S/0(y st ap )]£•(*) and E = N obj T. With these 
definitions Equation |^ can be written as, 
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N exp = E £{t)-,dl 
Jo at 



(11) 



A similar rc-dcfinition may be performed on Equa- 
tion [l] for the optical depth. The advantage of using 
E and £(t) instead of and £*{t) is twofold. Firstly, 
the exposure in object-years, E, is known accurately. 
Secondly, this substitution leaves £(t) containing the 
only reference to V^ t0 p (implicitly) and we can inves- 
tigate its convergence with magnitude easily. It is 
important to note that with the above definition of 
£ (t), the efficiency is no longer bound to lie below one. 
This slightly un-intuitive result is due to the fact that 
SyO(Kjtop) may be quite large (Table |l|), though in 
practice £(t) always lies below one. A practical way 
of viewing £ (t) is: given S/0(V s t op ) events, £(t) is 
the expected number of detected events for the given 
i. In the limit S/0(V stop ) ->■ 1.0 the efficiency £{t) 
recovers its usual meaning. 

Figure |ll| illustrates the integration of e(i,V), 
panel (b), over a LF, panel (a). The resulting function 
£{t) = [S/0(V < 24)] * £_(*) is shown in panel (c). 
Our two LFs, LF\ and LF2, and their corresponding 
values of S/0(V < 24) are displayed as dotted and 
solid lines, respectively. An important result is that 
£(t) is fairly robust to uncertainties in the LF fainter 
than V ~ 22. For durations less than 75 days the 
difference in £(t) as derived using either LF is less 
than 1%. The difference, however, becomes progres- 
sively larger for longer durations due to the relative 
difference in contributions from faint stars exhibited 
by the two LFs. At 300 days the difference is ~3%. 
In: " 
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we chose to favor LF2 over LF\ because of our 
HST completeness tests and because it seems unlikely 
that the LF will continue to rise so steeply beyond the 
clump for so long. Evidently any moderately differ- 
ent LF would produce only a small difference in the 
overall efficiency. 

In Figure fl2| we present the efficiency £(i) for selec- 
tion criteria set A and B using LF LF2. Also shown 
for comparison are the photometric efficiencies used 
in A96 and A97. Many of the differences between 
the old and new results, as seen in the sampling ef- 
ficiencies (§ 3.4), can also be seen here. Again the 



most striking difference is the much larger efficiency 
at long durations. As described in § 3.4 this has mul- 
tiple causes, including (1) almost three times more 
baseline in the A00 data than in the A97 data, (2) 
a longer duration cut of t < 600 days and (3) 6 high 
density fields which were previously split into two sep- 



arate years of data (due to an early generation of tem- 
plates) have been recombined. A new reason, unique 
to the photometric efficiency, is the contribution of 
faint stars to the efficiency for long duration events. 
This effect was troublesome in the A96 and A97 re- 
sults since we lacked adequate numbers of faint stars 
(only a few with V ~ 22 and none with V > 22) 
and thus were only confident that the photometric 
efficiency had converged for durations less than 150 
days. As we show below the A97 results were likely to 
have converged only for durations less than 100 days. 
However, this resulted in an underestimate in the ef- 
ficiency of less than 10% for durations around 150 
days. Criteria set A is closer in design to the A97 cuts 
and, as a consequence, follows the older results more 
closely for very short durations up to about 60 days. 
As noted above the sharp cut-off at durations around 
600 days seen in the sampling efficiency is smoothed 
over in the photometric efficiency because blending 
scatters intrinsically long duration events to shorter 
measured durations. 

5.2. Convergence 

Figure [l3| shows the convergence of £{t) with V s t op 
for the four combinations of selection criteria sets A 
and B and LF\ and LFi- In each panel the conver- 
gence of four different durations (t = 50, 100, 300, 
and 1000 days) is shown. The ordinate is in arbitrary 
relative units. Criteria set B is our loosest set of cuts 
and is the least convergent of the two selection crite- 
ria sets used in A00. An inspection of the figure also 
shows that LF\ gives a somewhat slower convergence 
than LF2. This is not surprising as LFi contributes 
a substantially larger number of fainter stars as com- 
pared with LF2. Criteria set A with LF2 converges 
for durations less than 1000 days by V^ top ~ 24.5 and 
criteria set B with LF2 converges for durations less 
than 300 days by V s t op ~ 24. If the true LF is closer 
in form to LF\ then it is possible that our efficiency 
results are somewhat underestimated for long dura- 
tion events (see § 5.3 below). However, since none of 



the 17 candidates in A00 have durations longer than 
300 days we feel that the convergence of the efficiency 
is more than adequate. 

5.3. Error Budget 

There are a number of potential sources of error in 
estimating our efficiency and we list the most impor- 
tant ones in Table ^. We crudely classify the errors as 
'signed' or 'unsigned'. That is, if we are reasonably 
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certain the effect would only increase (or decrease) 
the efficiency we classified it as 'signed' and marked 
it with the appropriate sign in the table. If we are 
not certain of it's sign we called it 'unsigned' and left 



a qrestion mark in the sign column of the table. 



The first unsigned error in the table (#1) is simply 
the uncertainty in our normalizations and was esti- 
mated to be ~ 20% in S 4.1.2. That is we are confi- 



dent that we know the number of stars to the limit 
V = 24 in our LMC fields to within - 20%. This is by 
far our largest source of error in the efficiency. The 
second unsigned error (#2) is an uncertainly in our 
incomplete knowledge of the shape of the LMC's LF 
at faint magnitudes. In 
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we used the difference 
between LF\ and LF 2 to estimate a 1-3% effect de- 
pending on the duration of the event. Another source 
of unsigned error (#3) is due to the finite number of 
Monte Carlo events, which we estimate to be < 1% 
based on the binomial distribution and the number of 
Monte Carlo events used. 

The fourth entry in Table || (#4) is signed and rep- 
resents a potential over-completeness bias. Because 
blending might over-populate objects with respect to 
stars in our calibration bin (17.5 < V < 18.5), this 
could cause an overestimation of £{i). In § 4.1.2| we 
estimated this could be as large as ~ 9%, but likely 
to be somewhat smaller (perhaps 6-8%). The next 
signed error (#5) is related to our choice of using 
only the nearest neighbor to each artificial star, when 
in fact as many as ~ 3% of PRFs had 2nd or 3rd 
nearest neighbors which were more sensitive to the 
lensed flux. Of course what matters is how many 2nd 
or 3rd nearest neighbors would have been detected 
but were not because the closest neighbor was used 
instead, and this is probably smaller than the number 
of 2nd and 3rd nearest neighbors who are more sen- 
sitive. The efficiency is likely to be underestimated 
somewhat by this error. The next two signed errors 
are small and relate to (#6) our use of a 'universal' 
LF (especially in the outer LMC disk fields where the 
true LF likely turns over at brighter magnitudes) and 
(#7) a possible slight underestimation of £ (t) if con- 
vergence has not been reached by V — 24 (e. g., if the 
the tr ue u nderlying luminosity function is closer to 

We leave two potential sources of error unexplored 
in this paper because of their highly uncertain and 
complicated nature. The first of these (#8 in the ta- 
ble) is the difficulty of incorporating binary star sys- 
tems into the determination of the LF. Binary stars 



could increase our exposure (thereby increasing £(t)) 
but also shift the LF toward fainter magnitudes (de- 
creasing £(t)). The effect on the LF of the local neigh- 
borhood is still controversial flKroupa 1995 , Reid & 



Gizis 1997) and any hope of resolving this for the 



LMC is beyond the scope of this paper, and left for 
future work. Another unexplored source of error (#9) 
is due to our lack of exotic lensing Monte Carlo events. 
The behavior of this source of error is difficult to es- 
timate as the shape of exotic lensing is likely to lower 
our effici ency, but caustic crossings might well in - 
crease it flDi Stefano fc Pcrna 1997j pi Stcfano 1999|) . 
Binary lensing events are likely to dominate this ef- 
fect, but the distributions of their event parameters is 
wide open to speculation (but see Alcock et al. 1999c). 
We plan to investigate the effects of a binary lens pop- 
ulation on our efficiency in a future paper. 

Given the generally small size of the individual 
signed errors relative to the size of the unsigned 
errors, as well as the Poisson-like 'counting' errors 
in the optical depth estimates (13-17 events gives 
~ 30 — 40% errors in the optical depth; see Table 10 
in A00), halo mass fraction, etc., and the fact that 
a number of these signed errors are of comparable 
size but opposite sign, we have chosen not to attempt 
to correct for them. Our estimates of their size are 
rough, but probably accurate to within a factor of 
two. There is no correct way to total systematic er- 
rors - but assuming they are uncorrelated we can sum 
them in quadrature to get a total of ~ 22%, not in- 
cluding errors #8 and #9, which may be the largest, 
but for which we do not have good estimates. Alter- 
natively, we can add the signed errors algebraically 
to get ~ 5%, then to this add the unsigned errors in 
quadrature, giving ~ 21% total error. So our best 
guess at the size of our error is 21-22% plus the un- 
known effects of binary corrections to the LF and ex- 
otic microlensing. 

5.4. Efficiency Differences Across Fields 

Our efficiency for each field is shown in Figure |l4|. 
These have been calculated in the manner described 
above and normalized with each field's S/O (Table 0). 
By splitting our ~ 1.2 million Monte Carlo events into 
30 fields, the shot noise has increased to 5-10%, which 
is somewhat field dependent. In principle the shot 
noise can be reduced by simply running more passes, 
but the CPU time required is somewhat prohibitive at 
the present time. Nevertheless the efficiency for each 
field is important when looking for gradients in opti- 
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cal depth or rates across the face of the LMC, which 
the present level of accuracy is adequate for, given the 
small number of events (13-17) seen toward the LMC. 
Differences seen between fields in Figure [Ti] are due 
primarily to (1) different normalizations for each field 
(Table [l]), physically corresponding to larger numbers 
of stars and thus probability of detecting events in 
some fields, (2) different sampling rates due to ob- 
serving strategy and weather across fields, and (3) to 
a smaller degree more crowded fields tend to have a 
larger effective area (area covered by stellar PSFs in 
the image plane) in which to detect the presence of 
magnified events, because of blending. 

5.5. Parameter Bias 

As discussed in the introduction, blending induces 
a bias in the measured parameters A max and t. We 
can quantify this bias using our Monte Carlo events, 
and have developed a method of statistically correct- 
ing for t bias when calculating optical depths. Fig- 
ure |l5| plots the ratio A' max /A max versus t'/i for a 

random samplp (N ~ 1 f) 5 ) of Monte Carlo pypnt.s- 



parlpl (a.) is for sampling pffipjpnpy pypnt.s and (h) is 



for photometric efficiency events (both use selection 
criteria, set R) The majority of sampling pfficipn 



events are rpmvprpH with little or no parameter bia.s 



in it (There is a. small amount of bia.s in (both 



directions) caused by sampling and weather gaps that 
effectively obfuscate very high magnification events.) 
As the sampling efficiency does not simulate blend- 
ing this result is not surprising. The results for the 
photometric efficiency events is quite different. A sig- 
nificant amount of blending can be seen in panel (b) 
of Figure [l5[ were approximately 40% of the Monte 
Carlo events are blended with an additional 10% or 
more of un-lensed flux. The trend of decreasing A max 
with decreasing t for highly blended events is well de- 
lineated in the figure. Note that the photometric effi- 
ciency events sample the CMD (luminosity function) 
in a square-root fashion (§ L2) and thus dim stars are 
underpopulated in panel (b) of Figure |l5[ Factoring 
this in, as is done below, further increases the amount 
of blending and the size of any potential correction. 

We can calculate the efficiency as a function of 
various parameters (e.g., u' min , u min , V, V obi , etc.) 
for the photometric efficiency, and some of these are 
shown in A00 (for example, see their Figure 6). How- 
ever, since the Monte Carlo events sample the CMD in 
a square-root fashion we must bin and re-weight with 
the correct LF. We must also bin and re- weight with 



a realistic distribution in t as the Monte Carlo events 
are added with durations uniform in log(t), which is 
unlikely to be the true distribution. We chose to re- 
weight with the distribution of durations predicted 
by a standard halo model with delta function mass 
0.5 M G ( priest 1991| ), since this distribution closely 
matches the data (average duration i ~ 92 days; see 
A00). Fortunately, model dependency is weak as the 
distributions in u m i n , weighted and un- weighted (that 
is using the log(t) distribution) arc quite similar. 

If the t bias is left uncorrected the optical depth 
will be underestimated (Equation ^) . There are two 
very different, but complimentary, ways to correct 
for this bias in T meas . Ideally one could perform mi- 
crolensing fits on individual events that allow for un- 
lcnscd flux in each passband. Then each event could 
be corrected for blending separately and the maxi- 
mum amount of information in the set of i's would be 
retained. In practice this approach is difficult as the 
lightcurves are not always well sampled and can have 
considerable photometric scatter making the blend 



fits uncertain flAlcock et al. 1996| , |Pratt 1997| , |Woz 



niak 1997, Han 1998). This situation can be greatly 



improved with accurate, dense, follow-up photometry 
on alerted events as is being done now routinely (Al- 

1997c 



brow et al. 1998a, Alcock et 



Rhie et al. 



1999). However, not all events are alerted. In addi- 
tion, blending fits are not always unique and only 
on fairly high magnification events can the fit pa- 
rameters be extracted with confidence. A second ap- 
proach is to correct the t bias in a statistical fashion. 
Since we have a priori knowledge of the distribution 
of source stars from the luminosity function <p(V) and 
can estimate the distribution of uncorrected durations 
alT/dt(t) we can compute a first order average correc- 
tion as, 



_ _ / dV J alt a(t, V) e(t, V) §(t) <P{V) 
JdV Jdie(t,V)§ (t)ct>(V) 



(12) 



where e(t, V) is defined as above and a(t, V) is defined 
as: a{t, V) — median(t' /i) or the median t bias for 
the Monte Carlo events in a bin of t and V. The 
quantity a(i, V) lies between zero and one and is a 
strong function of t and V in the sense that a(t, V) 
is one for bright events and approaches zero for both 
dim events and for long duration events. The results 
of Equation [l^ for criteria set A and B using LFs LF\ 
and LFi are tabulated in Table pi assuming mean 
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durations of 41, 92, and 130 days (corresponding to 
delta-function masses of 0.1, 0.5, and 1.0 Mq in the 
model dT/dt(t), respectively). The mean duration of 
the microlensing candidates in A00 corresponds most 
closely with our value of t ~ 92 days or an average 
mass of Mq ~ 0.5. This correction may be used 
in Equation [l] by simply substituting t' — > if /a to 
obtain, 

Note the optical depth defined in this manner does 
not scale simply as 1/cf. We have tried computing a 
as a function of measured parameters (such as Vobj 
and t') but find that it is fairly constant over the 
measured ranges. The use of a model distribution 
in t is somewhat worrisome. However, the range in 
corrections a for different t distributions is acceptably 
small, about 5% (Table |). The error induced in the 
optical depth is somewhat smaller (~3%; see Table 9 
in A00). 

To check this statistical correction we ran a set 
of Monte Carlos on a toy model using our artificial 
events. The Monte Carlos work as follows: (1) first 
we bin the recovered artificial events in t and V , (2) 
we create N exp 'events' by sampling the LF 4>(V) and 
t distribution function dT/dt for a simple halo model 
with a delta function mass 0.5 Mq (model S in A96, 
A97, and A00), (3) for each of these N exp 'events' we 
chose an artificial event to represent it by selecting 
at random from the appropriate bin in t and V, (4) 
three different optical depths for these N exp 'events' 
are computed using three different t's (the true t, the 
fit t' and the statistically corrected t st = t' /a). The 
number of N exp 'events' corresponds to the predicted 
number of observed events as computed by Equa- 
tion ^. The experiment is repeated 1000 times and 
the results for selection criteria set B and LF2 are 
displayed as histograms of r meas in Figure |l(| The 
dashed histogram corresponds to calculating T meas us- 
ing the fit t"s and clearly underestimates the optical 
depth {j mo dei = 4.7 x 10~ 7 for Model S; the ver- 
tical line in the figure). Using the true durations 
i to compute r meas gives the correct optical depth 
{solid histogram). The optical depth computed us- 
ing the statistically corrected durations, t' /a, is dis- 
played as the dotted histogram and correctly predicts 
the true optical depth in a fairly unbiased manner. 
Similar unbiased results are found using selection cri- 



teria set A, LFi, LF2, and using the dT/di distribu- 
tions with delta- function masses of 0.1 and 1.0 Mq. 
These Monte Carlos give us confidence that this sta- 
tistical correction can be used in an unbiased manner 
to compute the optical depth toward the LMC. 

6. SUMMARY 

The results of A00 for 5.7 years of photometric data 
toward the LMC rely critically on how well we un- 
derstand the detection efficiency. In this paper we 
calculated these efficiencies, correcting the most im- 
portant systematics effects with realistic models. To 
account for the wide range of issues due to blending 
we perform artificial star tests on a broad range of 
images. These artificial stars tests allowed us to em- 
pirically account for blending and to model our pho- 
tometry code's systematics. Correcting for blending 
also required an accurate knowledge of the LMC's LF 
because microlensing magnifies flux and increases the 
survey's sensitivity to dim, unresolved stars. 

We found that our sensitivity £(i) to dim, un- 
resolved stars vanishes for magnitudes fainter than 
V ~ 24 and durations less than i ~ 300 days. Our 
sensitivity in previous results (A96 and A97) was 
somewhat underestimated for durations greater than 
f ~ 100 days. We also found the t bias, due to 
blending, is of the order of 20% and we presented a 
method for statistically correcting this bias in our op- 
tical depth estimates. The method is complimentary 
to using blended t fits. 
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Fig. 1. — The 5.7-year sampling efficiency (normalized to u m in < 1-0) results for criteria set A (solid line) and 
criteria set B (dotted line). Also shown are the sampling efficiencies from A96 (long dash) and A97 (short dash) 
for comparison. 



24 



2.5 



Field 2 
Field 11 
Field 13 



.•7\ 



// 



Co 
+j 
CO 



o 
o 



0.5 



I/ ylN* 



if 

Nil 



: : I 



•\ 

'i !\ 



w V < 

A f: : ii 

Jj lJJ ! I 



16 



18 



20 



22 



24 



V 



Fig. 2. — Three HST LFs of MACHO fields are plotted with line codings displayed in the figure. Details of the 
relative normalizations can be found in the text. Note the good match in shape from V '~ 17 to T^~22. 
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Fig. 3. — Six Olsen (1999) HST LFs are plotted with line codings as displayed in the figure. Details of the relative 
normalizations can be found in the text. Note the good match in shape from V^17 to V^22. 
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Fig. 4.— Combined MACHO HST LF (dotted line) is plotted and compared with the combined Olsen (1999) HST 
LF (solid line). Also shown is the ground-based LF for MACHO field 13 (dashed line) which have been normalized 
as described in the text. The agreement in shape over the range 17 < V < 21 is good. 
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Fig. 5. — Four sample MACHO ground-based LFs (solid lines) ranging from low to high stellar density are labeled 
in the figure. LFi (histograms) has been offset and normalized (Table [|) to each field. The agreement in shapes 
between the individual ground-based LFs and LF2 is quite reasonable. 
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Fig. 7. — Four photometric response functions (PRFs) are plotted as recovered magnification A' versus input 
magnification A (open circles) and compared with the "sampling efficiency" analytic response function described 
in § H (solid lines). The dashed line is a least squares fit to the PRF data and in each panel the fit slope (fraction 
of 'lensed light' or blending fraction fuend) is labeled along with the distance (S in arcsecs) of the 'lensed' star to 
the nearest photometred object. 
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Fig. 8. — The four panels here correspond to the same PRFs shown in Figure |7|. The characteristic change in 
the relative error (error in magnitudes) A a is plotted versus A. Overplotted on each panel is the purely Poisson 
response assumed in the sampling efficiency (solid line; see § [|) and the behavior ~1/A (dashed line). See the text 
for more explanation. 



31 



< 

T3 



< 

X5 



0.5 





V ~ll 

' blend 1 S 



19\ 20 



21 



22 



60 



50 



40 



30 



20 



10 



> 



A = 1 .75 



16 



18 



20 



22 



24 



V. 



prf 



Fig. 9. — Plot of effective blend fraction (fit m = dA rec /dA) versus V magnitude for the <~ 59, 000 PRFs contained 
in the PRDB. Blending is not limited to faint stars, but clearly becomes worse for fainter magnitudes. Note that 
the distribution of PRFs shown here samples the LF in a square-root fashion (for reasons discussed in the text) 
and as such dim stars shown here are under-weighted. The family of curves and the right-hand axis are explained 
in the text. 
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Fig. 10. — Displayed here are four example Monte Carlo lightcurves with identical input event parameters (A max — 
2.2 and t = 80 days) generated using the same four PRFs displayed in Figures and ||[ The solid line is the standard 
microlensing fit described in § ||. The error bars are ± la. 
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Fig. 11. — Panel (a) displays luminosity functions LF\ {dotted line) and LF2 (solid line), which are used to integrate 
e(t, V) over V. The function e(t, V) is shown in panel (b) as a contour plot, with the contours corresponding to 
efficiencies of 0.001, 0.01, 0.1, 0.2, 0.3, 0.4, and 0.5. The function e(t,V) is described more fully in the text. The 
result of integrating e(t, V) over the two LFs and scaling by the appropriate normalization to obtain £ (t) is shown 
in panel (c). The line coding corresponds to that in panel (a). Note that £(t) is fairly insensitive to the shape of 
the luminosity function for dim stars because, (1) the function e(i, V) is small here and (2) although the faint stars 
are weighted less in LF 2 they also contribute less to the normalization S/0(V < 24). 
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Fig. 12. — Microlensing detection efficiency (normalized to u m i n < 1) for the 5.7-year MACHO data, as a function 
of event timescale t. The solid line sho ws t he 'photometric' efficiency computed for criteria set A, and the dotted 
line for criteria set B as described in 8 5.1 



For comparison the corresponding curves from A96 and A97 are also 



shown. 
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Fig. 13. — Convergence of £(i) with Vkop is shown for several different durations, t, and for the four combinations 
of selection criteria sets A and B and luminosity functions LF± and LF2 . 
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Fig. 14. — The photometric efficiency for each field is displayed for the 30 fields analyzed in A00. The different 
sampling rates, stellar densities, and crowding all contribute to the differences seen between fields. 
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Fig. 15. — Parameter biases in the fitting code are induced by blending as illustrated here in a scatter plot of 
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The table lists the normalizations (columns 7-10) for each of the 30 fields, along with conservative 
estimated errors. The columns 2-4 list the airmass, seeing (in arcsecs), and sky level (in electrons) for the 
MACHO template images. Column 5 is the offset for each field used in § 4.1.2j . Column 6 is the average 
density of each field in Objects /□' '. 
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Table 2 
PDRB Chunks 



ID 


Field 


Chunk 


Template 


Objs/W 


^* cond 


N PRF /grid 


NpRF 


1 


2 


9 


5672 


143.9 


67 


64 


3226 


2 


7 


40 


5668 


196.4 


66 


81 


4209 


3 


19 


22 


4562 


142.1 


72 


64 


3219 


4 


22 


38 


1696 


93.7 


64 


64 


3419 


5 


55 


54 


5311 


98.1 


69 


64 


3294 


6 


77 


27 


4203 


257.4 


67 


144 


6884 


7 


77 


50 


2962 


308.9 


66 


156 


7622 


8 


78 


6 


5663 


318.7 


69 


156 


8381 


9 


78 


47 


5663 


269.5 


77 


144 


8074 


10 


81 


18 


5673 


202.9 


72 


156 


6653 



The ten 5'x5' chunks used to create the PRDB span the observed range of density 
{Objs/W) and are widely spaced over the face of the LMC. For reference we give 
the template observation number, the number of observing conditions (N con d) used 
for each chunk, the number of artificial stars per grid (N pp F / 'grid) added to each 
chunk (60 grids per chunk) and the number of useful PRFs generated (Nppp) from 
each chunk. 
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Table 3 
Efficiency Table 



Bin 


tmin 


tmax 


„ 
t 


Criteria set A 


(~ -\ • J • „ „ i_ T") 

Criteria set B 










S(t) 


£{t) 


1 


1.00 


1.46 


1.20 


8.7e-04 


1.9e-04 


2 


1.46 


2.00 


1.70 


2.1e-03 


5.7e-04 


3 


2.00 


2.70 


2.32 


6 5e-03 


2 2p-03 


4 


2.70 


3.80 


3.20 


1.3e-02 


5.7e-03 


K 



o.ou 


^ 9D 


A AA 






6 


5.20 


7.10 


6.07 


5.2e-02 


3.7e-02 


7 


7.10 


9.70 


8.29 


8.7e-02 


7.8e-02 


8 


9.70 


13.30 


11.35 


l.le-01 


l.le-01 


9 


13.30 


18.20 


15.55 


1.5e-01 


1.8e-01 


10 


18.20 


24.90 


21.28 


1.9e-01 


2.5e-01 


11 


24.90 


34.20 


29.18 


2.5e-01 


3.2e-01 


12 


34.20 


46.80 


40.00 


2.8e-01 


3.8e-01 


13 


46.80 


64.20 


54.81 


3.2e-01 


4.4e-01 


14 


64.20 


88.00 


75.16 


3.6e-01 


4.7e-01 


15 


88.00 


121.00 


103.18 


3.9e-01 


5.1e-01 


16 


121.00 


165.00 


141.29 


4.1e-01 


5.4e-01 


17 


165.00 


227.00 


193.53 


4.2e-01 


5.3e-01 


18 


227.00 


311.00 


265.70 


4.3e-01 


5.3e-01 


19 


311.00 


426.00 


363.98 


4.1e-01 


4.9e-01 


20 


426.00 


584.00 


498.78 


3.2e-01 


4.1e-01 


21 


584.00 


800.00 


683.52 


1.0e-01 


1.7e-01 


22 


800.00 


1096.00 


936.37 


4.7e-02 


9.1e-02 


23 


1096.00 


1502.00 


1283.04 


3.3e-02 


6.8e-02 


24 


1502.00 


2000.00 


1733.14 


3.0e-02 


6.1e-02 



Efficiencies presented in A00 and also shown in Figure [12| for selection 
criteria set A and B. There are 24 bins in duration with t m i n and t max 
delineating each bin. The quantity t marks the logarithmic center of each 
bin. 
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Table 4 
Error Budget 



Number 


Sign 


Size 


Note 


1 


? 


-20% 


LF normalization/exposure 


2 


? 


-1-3% 


LF for faint stars V > 22.5 


3 


? 


<1% 


Finite number of Monte Carlo events 


4 


+ 


~6-8% 


Overcompleteness bias 


5 




-3% 


1st neighbor PRF bias 


6 


+ 


-1% 


LF in outer 'disk' fields 


7 




<1% 


Convergence for i > 100 days 


8 


? 


?% 


Binary stars in LF 


9 


? 


?% 


No exotic lensing events 



This table details our error budget for the efficiency. Column 
2 gives the probable direction, if known, in which the error might 
propagate to the efficiency, e. g., #4 could bias the efficiency to higher 
values, causing the optical depth to be slightly underestimated. Two 
possible corrections that are not estimated here are #8 and #9 (see 
the text). 



Table 5 
i Bias Corrections 



Selection 


LF 




Median a 




Criteria 




0.1 M© 


0.5 M© 


1.0 M© 


A 


1 


0.864 


0.826 


0.812 


A 


2 


0.867 


0.831 


0.818 


B 


1 


0.826 


0.777 


0.759 


B 


2 


0.828 


0.785 


0.770 



This table lists t bias corrections which may be used 
to statistically correct the unblended fit t's. This sta- 
tistical correction is a function of the selection criteria 
set used, but also sensitive to the average duration of 
the observed microlensing events, here delineated us- 
ing a standard Galactic halo model (model S) with 
delta-function masses of 0.1, 0.5, and 1.0 M©, corre- 
sponding to mean durations of 41, 92, and 130 days, 
respectively. We recommend using LF 2 and 0.5 M©. 
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